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ABSTRACT 

The initial mass function (IMF) of the first (Population III) stars and Population II (Pop 
II) stars is poorly known due to a lack of observations of the period between recombination and 
reionization. In simulations of the formation of the first stars, it has been shown that, due to the 
limited ability of metal-free primordial gas to cool, the IMF of the first stars is a few orders of 
magnitude more massive than the current IMF. The transition from a high-mass IMF of the first 
stars to a lower-mass current IMF is thus important to understand. To study the underlying 
physics of this transition, we performed several simulations using the cosmological hydrodynamic 
adaptive mesh refinement code Enzo for metallicities of 10~^, 10~^, 10~^, and 10~^ Zq. In 
our simulations, we include a star formation prescription that is derived from a metallicity 
dependent multi-phase interstellar medium (ISM) structure, an external UV radiation field, 
and a mechanical feedback algorithm. We also implement cosmic ray heating, photoelectric 
heating, and gas— dust heating/cooling, and follow the metal enrichment of the ISM. It is found 
that the interplay between metallicity and UV radiation leads to the coexistence of Pop III 
and Pop II star formation in non-zero metallicity {Z/Zq >10~^) gas. A cold (T <100 K) and 
dense (p>10~^^ g cm~^) gas phase is fragile to ambient UV radiation. In a metal-poor {Z/Zq 
<10~^) gas, the cold and dense gas phase does not form in the presence of a radiation field of 
Fo~10~^ — 10~^ erg cm~^ s~^. Therefore, metallicity by itself is not a good indicator of the Pop 
III— Pop II transition. Metal-rich (Z/Z0>1O~^) gas dynamically evolves two to three orders 
of magnitude faster than metal-poor gas {Z / Z q<1{)~^) . The simulations including supernova 
explosions show that pre-enrichment of the halo does not affect the mixing of metals. 



Subject headings: cosmology: theory — galaxies: high redshift — stars: formation 



1. INTRODUCTION 

In the last two decades, cosmological simula- 
tions have become an important tool for theo- 
reticians to simulate the structure formation in 
the universe from the primordial density fluctu- 
ations. The cosmological model that currently 
matches the observations of the cosmic microwave 
background (CMB) best is the so-called cold dark 
matter (CDM) model. According to hierarchi- 



cal structure formation ( Press fc Schechter 19741 ), 
small clumps are the first to collapse and form 
structures, and then through mergers and the 
accretion of matter they build up larger struc- 
tures. Cosmological simulations, based on the 
CDM model of hierarchical structure formation, 
predict that the first stars (30—300 Mq) have 
formed at redshifts z ^ 20-3 0, in dark matt er ha- 
los with masses of - 10^ (|Tegmark et aTlll997 : 
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Abel et alJl2QQ2l:lBrQmm et al 120021 : lYoshida et al 



2003HOVShea &: NQrmanl[2007^ 

In the current literature, the first stars are 
called Population III (Pop III) stars. The metal- 
licity of Pop III stars is so low that metal cooling 
does not have any effect on their formation. In 
the simulations of the formation of Pop III stars, 
it has been shown that, due to the limited ability 
of metal-free primordial gas to cool, the initial 
mass function (IMF) of Pop III stars is a few 
orders of n iagnitude more massive than the cur- 



rent IMF ([Abel et al.l l2QQ2|; lYoshida et all I2QQ6I : 



Bromm et al. More recently, it has been 



suggested that Pop III stars are not necessarily 
that massive and can be in the range 10—100 Mq 



(iTumHnsonl I2OO6I l2007a||bD- 

On the other hand, 

from observations of the nearby universe we know 
that the pres ent d ay stellar mass scale is ^ 0.3 



Mq (|KrouDall2002l : IChabrierll2003h . 



The chemistry of zero metallicity (primor- 



thors (iDakarno & LcddI 198?!: Abe] 


et al.l 


1997; 


GaUi k Palla"l998': lAbel et al. 2001: 


Glover 


2008; 


Glover & Abel 2008:lTurk et al.l2011ah. Once the 



gas has virialized in the potential wells of dark 
matter halos, additional cooling is required for 
the further collapse of the gas and to form stars. 
The formation of a star depends on the ability of 
interstellar gas to cool and form dense molecular 
clouds. The modest cooling ability of primordial 
gas leads to high masses for Pop III stars. The 
cooling efficiency of star forming gas will be sig- 
nificantly affected by the addition of radiation 
and metals after the formation of the first stars. 
Hence, the chemical composition and the radia- 
tion environment of interstellar gas are the key 
parameters to study. 

The main coolants for primordial gas with tem- 
peratures T > 10^ K are Lya emission of neutral 
atomic hydrogen (H I, 1216 A), and ionized he- 
lium (He II, 304 A). Below this temperature, the 
dominant coolan t in primordial gas is m olecular 
hydrogen (H2). ISaslaw fc Zipov realized 



the importance of gas phase H2 formation in a 
primordial gas for the formation of protogalactic 
objects. At low densities (n < 10^ cm~^), H2 
can form via intermediate ions H2^ and H~. The 
cooling rate of H2 (Ahs) scales with density as 
Ah2 n^ at low densities (n < lO^cm"^), where 
radiative de-excitation dominates, and as Ahs cx n 



at high densities (n > 10^ cm~^), where colli- 
sional de-excitation dominates. Molecular hydro- 
gen does not have a permanent electric dipole mo- 
ment while the hydrogen deuteride (HD) molecule 
does, which makes HD a better coolant than H2. 
Although the primordial deuterium abundance is 
small relative to hydrogen (no/^H = 4xl0~^), 
chemical fractionation leads to an enhancemen t 
of the ratio, n nn/nn ^ ^ 10~^ (IPuv et al.l (Il993h : 
Galh fc Fallal (|l998l ): lstancil et al.l (|l998[ )). If pri- 



mordial gas is significantly ionized then HD cool- 
ing can lower the te mperature to the level of the 
CMB at z - 10-20 (| Johnson fc Bromml[2006l ). 

When the first stars form, the universe becomes 
much more complex. According to their initial 
masses, they will either explode as core-collapse 
supernovae (10 Mq < < 140 Mq) or as pair in- 
stab ility supernovae (PISN; 1 40 Mq < < 260 
Mq , iHeger fc WooslevI (|2002[ ) ) . Through these su- 
pernova explosions (SNe), the interstellar medium 
(ISM), and the intergalactic medium (IGM) will 
be enriched with metals. From that point on one 
has to take into account the cooling from the fine- 
structure lines of metals and the rotational tran- 
sitions of molecules like CO. Moreover, gas— dust 
heating/cooling can be important as well, in the 
sense that dust grains will be heated to at least 
the CMB temperature. Dust allows the efficient 
formation of H2 and HD (jCazaux fc Spaansll2004l . 
l2009h . and may collisionally heat / cool the gas de- 
pending on the si^ n of gas temperatur e minus 
dust temperature (^Schneider et al.l |2002[ ). Also, 
dust grains attenuate UV radiation. For z > 10, 
CMB photons can be an excitation source (ra- 
diative pumping) o f atomic and molecular levels 
([Smith et al.ll2009h . This allows the CMB tem- 
perature to act as a thermodynamic floor, below 
which gas cannot cool, provided that collisional 
de-excitation domin ates the remova l of populatio n 
from excited states (jSpaans fc Silkl (|2000l . l2005l )). 
Furthermore, SNe will initiate shock waves that 
propagate through the ISM. These shock waves 
can heat up the ISM and cause a delay in the for- 
mation of the next generation of stars or can com- 
press the gas, which makes it collap se and hence 



give ri se to further star formation dBromm et al 



2008D 



2OO3I: ISpringel fc HernquistI 120031 : IWhalen et al 



The aim of this work is to (1) compute at what 
metallicity a cold and dense gas phase emerges and 
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(2) assess the sensitivity of this phase to back- 
ground UV radiation. This paper is structured 
as fohows. In Section [21 we detail our cosmo- 
logical/hydrodynamic simulation. Specifically, we 
discuss the relevant cooling and heating processes 
for star formation and the global characteristics 
of early star formation in primordial galaxies. In 
Section [H we highlight the implications of our 
findings. Finally, in Section |4] we discuss our re- 
sults and present our conclusions. 



2. SIMULATIONS 

In this work, we use the cos niological adaptive 
mesh refineme nt code E nzo ("Brv an fc Normanl 
19971 : [O^Shea etld] l2QQ4l ). We perform simula- 
tions in a three-dimensional periodic box with a 
side length of 1 h~^Mpc, initialized at z = 99. The 
size of the root grid is 128^ with three nested sub- 
grids, each refined by a factor of two. The finest 
grid has an effective resolution of 1024^ with a side 
length of 125 h~^kpc. This resolution results in a 
dark matter and baryonic matter resolution of 2 
and 0.4 M©, respectively. Refinement is restricted 
to the finest grid and occurs during the simulations 
whenever the baryonic matter, or dark matter den- 
sity, is greater than the mean density by a factor 
of four or eight, respectively. The maximum level 
of refinement that is reached in the finest grid 
is eight. Refinement occurs such that the Jean 
length is always resolved by at least eight cells, 
this ensures that we meet the Truelove criterion, 
which requires the Jeans length to be resolved 



For the analysis of our cosmological simulations 
we use YT, a cross-platfo rm analysis toolkit writ- 
ten in Pvthon (|Turk.2QQ8jTurk et al.ll2Qllb[ ). 



by at least four cells on each axis (iTruelove et al. 
19971 ). The virial mass of our progenitor halo at 
redshift z = 21 is Myir 3.4 x 10^ M©, where Myir 
is the mass in a sphere that encloses an average 
dark matter overdensity of 200. We use Wilkinson 
Microwave Anisotro py Probe five-year cosmologi- 
cal parameters (Komatsu et al. 2009), which have 
the following values: = 0.7208, Qrn = 0.233, 
= 0.0462, as = 0.9, and h = 0.701. Here, Oa 
is the vacuum energy, Qrn is the matter density, 
fti) is the baryon density, as is the variance of 
random mass fiuctuations in a sphere of radius 8 
Mpc, and h is the Hubble parameter in units 
of 100 km s~^ Mpc~^. We focus on a single halo 
with a dark matter mass of ~ 10^ Mq at z ~ 5 
which is expected to be a typical mass for the halo 
population at that redshift. 



2.1. Cooling and Heating Processes in the 
ISM 

Metal enriched gas cools more efficiently by 
fine-structure lines of [C II] (157.74 /im), [O I] 
(63.18 /im, 145.5 /im), [Si II] (34.8 /im), [Fe II] 
(25.99 /im, 35.35 /im), and ro tational lines of 
CO t han by HD or H2 emission (jSantoro fc Shulll 
|2006[ ). In the outermost layers of a star- forming 
cloud, the so-called photon dominated region 
(FDR), temperatures can increase up to 1000 K 
due to a strong UV radiation field, and cooling 
results from the fine-structure lines of [C II] and 
[O I]. Deeper into the cloud, the temperature de- 
creases to ~ 30 K through the balance between 
cosmic ray and dust heating and low- J CO rota- 
tional line cooling. 

In our simulations, the gas is homogeneously 
pre-enriched to some non-zero metallicity at red- 
shift of z = 30. We perform a series of simulations 
using cooling models for four different metallic- 
ities, 10-"^, 10-^, 10-^ an d 10~^ Zc:,, derived 
from the chemical network of iMeijerink fc SpaansI 
(I2OO5I ). In their FDR code, they include cooling 
from fine-structure lines of carbon (C+ and C), 
oxygen (O), molecular lines from species like car- 
bon monoxide (CO), H2, HD, and water (H2O). 
All level populations are computed under statis- 
tical equilibrium and the chemistry includes gas 
phase and ^rain surface formation of H2 and HD 
( Cazaux fc SpaansI l2004i 12009 ) and line trapping 
usi ng: the multi-zone e scape probability method 
of iFoelman fc SpaansI (2005). The cooling ta- 
bles that we use, depend on FUV radiation field 
strength, ambient gas velocity dispersion, temper- 
ature, metallicity, and H2 abundance. As such, 
they provide an accurate treatment of the ther- 
mal and chemical balance of low- metallicity FDRs. 
The radiation field strength (Go) here is the fiux 
between 6— 100 eV and it extends beyond 13.6 eV 
according to a starburst99 spectrum for a Salpeter 
IMF, i.e., the H II region is computed as part of 
the FDR. Typically, the UV fiux diminishes very 
quickly beyond 13.6 eV. The cooling tables enjoy 
a range in irradiation (Go=10~^, 10~^, 10^, 10^, 
10^ 10^ 10"^), H2 abundance {/h^^ 10"^ lO""^, 
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10-^, 10-^ 10-\ 0.2, 0.4, 0.5) and metallicity 
(10-^ 10-^ 10-^ 10-\ l Z©). we adopt Milky 
Way like abu ndance ratios b ased on the values of 
lAsplund et al.l (j2005l ) and ljenkinsi (|2004l ). A sum- 
mary of our simulations is listed in Table [TJ 

At high column densities cooling and heating 
can be suppressed due to optical depth and dust 
opacity effects. In the case of an optically thick 
medium, a large part of the emitted photons are 
thus reabsorbed. This results in a lower critical 
density where the critical density is defined as the 
density at which the radiative de-excitation rate 
equals the collisional de-excitation rate. At den- 
sities much below the critical density (n <C ncr), 
radiative de-excitation dominates over collisional 
de-excitation, whereas at high densities (n > ncr), 
collisions dominate the de-excitation process and 
the gas is in local thermodynamic equilibrium 
(LTE). Radiation is trapped if the opacity is con- 
centrated on small physical scales which is the case 
in our simulations for the warm and dense gas. 
Since we do not resolve these PDR sub-structures 
in our simulations we approximate these opacity 
effects as follows. 

We adopt a turbulent coherence length L of 
0.3 pc and a local turbulent velocity dispersion 
of dV = 3 km s~^ to mimic the properties of a 
turbulent region. This yields a formal velocity 
gradient dV/L for scales larger than L. When 
the cooling is either optically thin (r <C 1) or 
optically thick (r ^ 1) the choice of dV and L 
is not important. For the radiative transfer in 
the cooling lines (under statistical equilibrium), 
we use the multi-zone esc ape probability code of 
Poelman fc SpaansI (|2005[ ) . Each zone in a PDR 



is treated this way, and a correction to the cooling 
and heating rates due to line t rapping and optica l 
depth effects is derived. See Wada et al. (|2009l ) 
for more details. 

For an incident UV radiation field, photo- 
electric emission from (small) dust grains and 
polycyclic aromatic hydrocarbons is the dominant 
heating source in the neutral ISM. Dust grains can 
absorb an FUV photon which leads to the ejec- 
tion of an electron, carrying some of the photon 
energy away in the form of kinetic energy. Then, 
through elastic collisions this excess kinetic en- 
ergy heats up the gas. In our simulations we add 
photoelectric heating through 



Here, Gq ^^^^ is the radiation field attenuated by 
dust absorption and is given by 



G, 



0,dust 



= Goexp(-1.8A,) 



(2) 



where Go 



a fiux of 1.6x10 ^ erg cm ^s ^ 



1, in Habing units, corresponds to 
Aj^ is the line- 
of-sight visual extinction at optical wavelengths 
caused by interstellar dust and the total num- 
ber density of hy drogen i s, nn = n(H) + 2n(}io ) 
([Bakes Tielens 1994; M eiierink Spaansll2005l ). 
e is the heating efficiency, dependent on Go , the ki- 
netic gas temperature , and the electron density 

1/2 

Ue as GqT^.' /ue, and is given by 



e = 



4.87 X 10-2 

[1 + 4 X 10-3(^!^)0.73] 

^ 3.65 X 10-2(^)°-^ 
[l + 2xl0-4(^)]' 



(3) 



On the other hand, gas and dust tempera- 
tures are not equal in the ISM. At high num- 
ber densities, uh >10 ^ gas— dust heat- 
ing/cooling becomes important and therefore we 
add gas— dus t heating /cooling into our simulations 
as given by (iMeiierink SpaansI I2OO5I ): 



Icon = 1.2 X 10 nH(-— -j ' ( ) ' 

— 75 

x[l-0.8exp( — )](rrf-rfc), 

J-k 



(4) 



where Td is the dust temperature, Tk is the gas 
temperature, and amin is the minimum grai n size 
which we take as 1 fiva (Nozawa et al."'2003). 

Deeper into a cloud the UV radiation is attenu- 
ated by dust and cosmic ray heating becomes im- 
portant. Ionization by cosmic rays pr oduce s en- 
ergetic electrons. Gla ssgol j^j^ Lan gerl (|l973l ) and 
Cravens fc Dalgarnol (|l978[ ) calculated that about 
8 eV of heat is deposited in a molecular gas per 
primary ionization. Includi ng helium ionization, 
Tielens fc HollenbachI (|l985[ ) find for the heatine: 
rate 



FcR = 1.5 X 10-^^Cn(H2) ergcm"^ s" 



(5) 



10 



-24 r^' — S - 

eGcdustHH erg cm s 



(1) 



where ^ is the cosmic ray ionization rate per H2 
molecule. In our simulations, we scale ^ with the 
derived star formation rate (SFR) of our simula- 
tions such that for an SFR of 1 Mq yr~^ and a 
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Table 1 
Simulation parameters 



Run 


L 


Metallicity 


SFFB 


Go 




(Mpc) 


{Zq) 




(Habing) 


Zl-Gl 


1 


10-^ 


No 


10~2 


Z2-G1 


1 




No 




Z3-G1 


1 


10~^ 


No 


10~^ 


Z4-G1 


1 


10-^ 


No 


10~2 


Zl-GlO 


1 


10-^ 


No 


10-1 


Z2-G10 


1 


10~2 


No 


10-1 


Z3-G10 


1 


10~^ 


No 


10-1 


Z4-G10 


1 


10~^ 


No 


10-1 


ZS2 


1 


10~2 


Yes 


10-2 


ZS4 


1 


10-^ 


Yes 


10-2 



Note. — Column 1: Simulation name; Column 2: 
Simulation box size; Column 3: Simulation pre- 
enrichment; Column 4: Star formation and feedback; 
Column 5: UV background strength. 



Salpeter-like IMF between 1 and 100 M0, we take 
C to be SxlO-^V^ (fSpaans Silk.20051 



2.2. Star Formation 

In order to understand the underlying physics 
of the transition from a high-mass IMF of Pop III 
stars to a lower mass current IMF we have run 
simulations for different metallicities and different 
UV radiation backgrounds and have incorporated 
the metal yields of SNe. 

In our simulations where star formation is in- 
cluded, a star particle is formed when a grid 
cell has a density exceeding 0.01 cm~^ and a 
temperature below 10^ K. We follow the hydro- 
dynamic transport of metals produced by the 
SNe to the enrichment of halos. We do not use 
the chemical network of Enzo, except for the 
non-equilib rium formation of H 2 (rate equation 



method of lAnninos et al. (|l997l )) but our A(Go5 



fH2^ Z/Z0) cooling tables instead. 

In order to compute the metal yield by the first 
stars we need to decide on some parameters like 
the mass of the halo, IMF of the first stars, star 
formation efficiency, and the number of SNe. From 
their simulations of pre-gala ctic structure forma- 
tion, Machacek et al.l (2001) determined that the 
minimum mass of a halo that hosts a massive pri- 



mordial star is 



exp(AL)(1.25xlO^ 
+8.7 X IO^fO^^i), 



(6) 



where Mmin is the minimum halo mass that con- 
tains a cold dense core; fed is the fraction of 
gas that is cold (T < 0.5 T^ir) and dense (p 
> 330 cm~^);,and Flw is the fiux within the 
Lyman— Werner (LW) bands in units of 10~^^ erg 
s~^ cm~^ IIz~^. We choose fed = 0-02 for a typ- 
ical UV background of J = 10~^^ erg s~^ cm~^ 
Hz~^ sr~^, so that a halo with a minimum mass 
of 4.16x10^ Mq can support cold dense gas which 
will eventually form a primordial star. Note here 
that J = J_2i corresponds to a total UV flux of 
Fo ^ 4x10-^ erg s"^ cm-^ or Go ^ 2.5x10-^. 
For a Salpeter like IMF, Go — 1 corresponds to 
an SFR of about 3 yr-^ (10 kpc)-^ 

High redshift, low mass galaxies are conjectured 
to have properties similar to those of local dwarf 
galaxies. In these galaxies the star f ormation effi- 



cienc y, f^, ranges from 0.02 to 0.08 (|Tavlor et al. 
1999[ ). In this work we take = 0.05. 



3. RESULTS 

3.1. Multi-phase ISM 

In our cosmological simulations for the Z2 — Gl 
and Z4— Gl runs we find (as expected) that gas 
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Density (g cm"') Density (g cm"') 

Fig. 1. — Density— temperature profile of the central 50 kpc of a halo at redshift 2; ~ 5 for metallicities of 
Z/Z0=1O-'^ (left) and Z/Z0=lO-2 (right) for a radiation field of Go=10-^ 




Density (g cm"') Density (g cm"') 

Fig. 2. — Density— temperature profile of the central 50 kpc of a halo at redshift z ^ 20 for metallicities of 
Z/Z0=1O-'^ (left) and Z/Zq=10-^ (right) for a radiation field of Go=lO-^. 
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in the high- met allicity case cools down to lower 
temperatures and has higher densities than in the 
low-metallicity case. This is shown in Figure [T] 
where we plot the density— temperature profile of 
a halo at redshift z ~5 for metallicities of Z/Z© = 
10-"^ (left) and Z/Z0=lO-2 (right) and a radia- 
tion field of Go=10~^. In the high-metallicity case 
there is more gas in the central 50 kpc. This is be- 
cause the cooling is more efficient, so that the gas 
has less pressure support against gravity and the 
halo accretes more material. When we look at the 
high densities (p >10~^^ g cm~^) in these plots 
it is clear that in the high-metallicity case there 
is a lot more gas at low temperatures than in the 
low-metallicity case. This gas is identified as the 
cold and dense gas phase of the multi-phase ISM. 
The key feature here is that the redshift at which a 
multi-phase ISM is established depends on metal- 
licity. This is shown in Figure [2] where we plot 
the density— temperature profile of the same halo 
as shown in Figure [1] at z ~20 for metallicities of 
Z/Zo = 10-"^ (left) and Z/Z© = IQ-^ (right). We 
see here that the cold dense phase is completely 
lacking in the Z4— Gl run whereas it is already 
present in the Z2— Gl run. This means that metal- 
rich gas cools more efficiently and therefore the 
multi-phase ISM is established at earlier times in 
the higher metallicity case than in the lower metal- 
li city case. These results are consistent with those 
of ljappsen et all ( 2QQ9h . although they did not in- 
clude background radiation. 

Due to the shorter cooling time a halo in the 
metal-rich case evolves dynamically faster and be- 
comes more compact at z = 5 compared to a 
metal-poorer halo. The effect of metallicity on 
a halo can be seen in Figure [3l which shows 
slices of the density, temperature, and the Jeans 
mass (Mj) of the central 1 kpc of a halo for the 
Z4-G1 (top), Z2-G1 (middle), and Zl-Gl (bot- 
tom) runs in the x — z (left column) and x — y 
(right column) planes at redshift 5. We compute 
the Jeans mass as 

M, = i^-^G^m^f/H^/^', (7) 

where p is the mass density, G is the gravitational 
constant, is the Boltzmann constant, T is the 
temperature, ja is the mass per particle, and ttih is 
the hydrogen mass. Because the high-metallicity 
gas {Z/Zq > 10~^) cools to lower temperatures, 
the Jeans mass in those runs is two orders of mag- 



nitude smaller than in the low-metallicity {Z/Zq 
= 10~^) case. We also see that the Jeans masses 
for metallicities Z/Zq = 10~^ and Z/Zq = 10~^ 
are comparable. This is because above a metallic- 
ity of about Z/Zq = 10 ~^ the cooling efficiency 
of ambient gas no l onger increases stron gly with a 
rise in metallicity ( Spaans fc Meijerinkll2QQ8[ ). In 
Figure HI we plot the density— temperature pro- 
file of the central 50 kpc of a halo at z ~ 5 for 
the Zl— Gl run. The phase diagrams of the runs 
Zl-Gl (Figure m) and Z2-G1 (Figure [B right) 
are indeed similar. 



3.2. UV Background Radiation 

In order to see the effect of a UV radiation field 
on the evolution of gas we performed multiple sim- 
ulations with different background radiation fields 
for the same metallicity. In Figure [5l we plot the 
density— temperature profile of the central 2 kpc 
of a halo at z ^ 8 for metallicities of Z/Zq = 10~^ 
(top) and Z/Zq = 10~^ (bottom) with two differ- 
ent background radiation fields of Go = 10~^ (left) 
and Go = 10~^ (right). When we increase the ra- 
diation field for a metallicity of Z/Zq = 10~^, 
the cold and dense phase does not form anymore, 
which is clearly seen in the phase diagrams of Fig- 
ure [5]: At densities above p = 10~^^ g cm~^ there 
is no cold gas present anymore in the Z3-G10 run 
(top right). This is because gas cannot cool down 
as efficiently as in the lower radiation field case 
due to more heat input and dissociation of H2 
and CO. Thus, the halo takes more time to col- 
lapse and become condensed. In our simulations 
we followed the halo until redshift z = 5, and in 
the Z3— GIO run the gas has still not collapsed at 
this redshift. On the other hand, in the Z/Zq = 
10 ~^ case the cold dense phase survives regardless 
of irradiation. This indicates that the cold dense 
phase is fragile to UV radiation for metallicities 
of Z/Zq < 10~^ and robust to UV radiation for 
metallicities of Z/Zq > 10"^ 

As is shown in Figure [6l where we plot the 
density profile of the central 2 kpc of a halo 
for the Z3-G1 (top left), Z3-G10 (top right), 
Z2-G1 (bottom left), and Z2-G10 (bottom right) 
runs, the effect of UV irradiation on the model 
galaxy morphology is large for Z/Zq = 10~^ and 
quite modest for Z/Z© = 10~^. In fact, we see 
that in the case of low metallicity and strong ra- 
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Z/Zo=10-^ 






Z/Ze=10-^ 




1 












i 


Z/Z0=1O-^ 

i 






Z/Zo=10-^ 




W^mk 


^^^^^^^ ' '^^^1 

^^w- ^^^^ 






Z/Z0=1O> 














^^^^ 

le-23 le-22 le-21 le-201e-19 
Density (g cm"^') 


10 100 1000 le4 
Temperature (K) 


1 10 100 1000 le4 le5 le6 
Jeans Mass (M©) 


le-23 le-22 le-21 le-20 le-19 
Density (g cm"^) 


10 100 1000 le4 
Temperature (K) 


1 10 100 1000 le4 le5 le6 
Jeans Mass (Mq) 



Fig. 3. — From left to right: Slices of density, temperature and Jeans mass of the central 1 kpc of a halo at 
redshift 5 for metallicities of Z/Z© = IQ-"^ (top), Z/Z© = IQ-^ (middle), and Z/Z© = 10"^ (bottom) in the 
x-z (left column) and x-y (right column) plane for a radiation field of Go = 10~^. 
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Fig. 5. 



- Density— temperature profile of the central 2 kpc of a halo at redshift 7.88 for a metallicity of 
10~^ with radiation fields of Go = 10~^ (top left) and Go = 10~^ (top right), and a metallicity of 
10~^ with radiation fields of Go = 10~^ (bottom left) and Go = 10~^ (bottom right). 
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Fig. 6. — Density slices of the central 2 kpc of a halo at redshift z = 7.88 for a metallicity of Z/Zq 
with radiation fields of Gq = 10~^ (top left) and Go = 10~^ (top right), and a metallicity of Z/Zq 
with radiation fields of Go = 10~^ (bottom left) and Go = 10~^ (bottom right) in the x — z plane. 
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diation field (Z3— GIO, top right) tiie compact 
disk-like structure is no longer formed, while in 
the higher metallicity and strong radiation field 
(Z2— GIO, bottom right) the disk-like structure is 
still present. This is because the high- metallicity 
gas is able to cool faster so that it collapses and 
reaches higher densities earlier in its evolution, 
and thus builds up column density to protect the 
cold and dense phase. 

In Figure [71 we plot the redshift at which a 
multi-phase ISM is established depending on the 
metallicity and the radiation field. We find that, 
to first order, the cold and dense phase survives 
if the ratio Fq/ Z is smaller than ~10~^ erg cm~^ 
s~^ in solar units. 

All this has a significant influence on the sub- 
sequent star formation. The Jeans mass sets a 
typical value for a cloud to collapse gravitation- 
ally and fragment, depending on both temperature 
and density of the medium. Although the metal- 
enriched gas has a higher cooling efficiency than 
metal-poor gas, under the influence of a back- 
ground radiation field the gas will stay hot and 
the Jeans mass will remain high. Therefore, the 
UV radiation background is a key parameter to 
take into account when studying the transition 
from Pop III to Pop II stars. In our simulations 
we see that including a constant radiation back- 
ground raises the critical density value for this 
transition from Zcr ^ 10~^-^Zq to Zcr ^ 10~'^Zq 
if Fo>10~^ erg s~^ cm~^. 



3.3. Dynamics 

When we look at the density slices of the 
Zl— Gl and Z2— Gl runs in Figure [3] in the x — z 
plane we see an elliptical spiral structure and in 
the X — y plane a flattened structure. This is a 
disk-like distribution in three dimensions and is 
more extended and massive in the Zl— Gl run 
(M = 5x10^ Mq) than in the Z2-G1 run (M = 
3.4x10^ ^o)- Because the gas cools faster in the 
halo of the Zl— Gl run, more gas is accreted and 
it settles faster into a disk than in the Z2— Gl run. 

In Figure [HI which is a volume rendering plot, 
we show the three dimensions density snapshots of 
the Zl-Gl (top), Z2-G1 (middle), and Z2-G10 
(bottom) runs for redshifts z = 7.88 (left), z = 
6.24 (middle), and z = 5 (right). We see that in 
the Zl— Gl run more structures are formed and 
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Fig. 4. — Density-temperature profile of the cen- 
tral 50 kpc of a halo at redshift z~5 for a 
metallicity of Z/Z0=1O~^ and a radiation field of 
Go=10-2. 
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Redshift 

Fig. 7. — Radiation field strength vs. redshift. 
Symbols mark the redshifts where a cold dense 
gas phase is established. The arrow represents the 
runs with Z/Zq=10-^ and Z/Zq=10-^, for a ra- 
diation field of Go=10~^, which do not develop a 
cold dense gas phase before z = 5. 
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Fig. 8. — Three-dimensional, volume rendered, density snapshots of the runs for metallicities of Z/Zq = 
10~^ (top) and Z/Zq = 10~^ (middle) for a radiation field of Go = 10~^ and Z/Zq = 10~^ (bottom) for a 
radiation field of Go = 10~^ at redshifts z = 7.88 (left), z = 6.24 (middle), and z = 5 (right). Size of the 
box for redshifts 7.88 and 6.24 is 60 kpc, and for redshift 5 is 10 kpc. 
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generally they have higher densities than in the 
Z2-G1 and Z2-G10 runs. When we look at the 
middle and bottom plots we clearly see the effect 
of UV irradiation on the evolution of a halo. In 
the Z2— GIO run, structures are more fluffy and 
take more time to become dense and compact 
due to the heat input. In these plots fluffiness is 
represented by the spatial extent of the different 
colors. Although the halo in the Zl — Gl run expe- 
riences multiple and more violent mergers as well 
as a higher rate of gas accretion, due to the faster 
cooling the halo at redshift 5 is denser and more 
compact than the Z2 — Gl case. In Figure [9l we 
plot the amount of mass for the densities indicated 
along the x-axis and for redshifts z = 7.88 (left), 
z=6.24 (middle), and z = 5 (right). In these plots 
it is clearly seen that there is more mass at higher 
densities in the Zl— Gl run than the Z2— Gl run. 

In Figure [TOl we plot the fraction of the halo gas 
mass that has a Jeans mass (Mj) below a certain 
threshold value for the Zl — Gl (black— dashes) 
and Z2— Gl (red— dots) runs. The profile of the 
low-metallicity case starts with a much steeper 
rise at low Jeans mass thresholds but then levels 
off for higher Jeans mass thresholds more rapidly 
than the high metallicity case. This initial steep 
rise seems to be due to a lack of low Jeans masses 
(< 30 Mj) in the central 1 kpc halo of the Z2-G1 
run. The fact that at any given threshold there is 
always a larger fraction of gas in the Zl — Gl run 
shows that the Jeans mass is generally lower in the 
high metallicity case than in the low-metallicity 
case. 

As previously mentioned. Figure [3] shows a 
disk-like structure for metallicities Z/Zq = 10~^ 
and Z/Zq = 10~^. In order to see if this disk- 
like structure is stable we plot the time evolution 
of specific angular momentum (see Figure [TT]) . 
In this plot we see that the specific angular mo- 
mentum increases with time indicating that the 
disk-like structure becomes more rotationally sup- 
ported. We also see several peaks which are caused 
by recent merger events. 

Next, we plot the spin parameter evolution 
of our main halo for the Z4— Gl run (see Figure 
[12]) . The spin parameter represents the degree 
of rotational support available in a gravitational 
system. We compute the dimensionless spin pa- 
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Fig. 9. — Gas mass vs. density threshold of the 
central 10 kpc for metallicities of Z/Zq = 10~^ 
and Z/Zq = 10~^ for a radiation field of Go=10~^ 
and Z/Zq = 10~^ for a radiation field of Go = 
10-^ at redshifts z = 7.88 (top), z = 6.24 (middle), 
and z = b (bottom). 
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Fig. 10. — Gas mass fraction vs. Jeans mass 
threshold for metalhcities of Z/Zq = 10~^ (black 
dots) and Z/Zq — 10~^ (red dashes) at redshift 
z = 5 for the central 1 kpc. 
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Fig. 11. — Specific angular momentum vs. red- 
shift for metalhcities of Z/Zq — 10~^ (black dots) 
and Z/Zq = 10~^ (red dashes) for the central 1 
kpc. 
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where L is the angular momentum, E is the en- 
ergy, M is the mass of t he object, and G is 



Barnes fc Efstathiou 



the gravitational constant. 
(1987) have calculated that the mean spin param- 
eter is < A > = 0.05 and it has been shown that 
halos which have suffered a recent major merger 
tend to have a highe r spin paramet er than the av- 
erage (e.g., Hetznec ker fc BurkertI (|2QQ 6)). As is 
shown in Figure [121 the spin parameter of our halo 
peaks around redshift z = 24, where a recent ma- 
jor merger occurred, and overall it has a value of 
^ 0.04. 

In Figure [131 we show the mass radial pro- 
files of the Zl— Gl and Z2— Gl runs, top and bot- 
tom, respectively. In both runs, the inner parts 
of the halo are dominated by baryonic matter. 
On the other hand, in the Z2 — Gl run dark mat- 
ter takes over at a smaller radius (700 pc) than 
in the Zl— Gl run (950 pc). This is because in 
the Zl— Gl run the halo is larger and more mas- 
sive baryonically than in the Z2 — Gl run. Once 
more, this shows that the halo evolves dynami- 
cally slower in the low metallicity run. 

In Figure [HI we plot the mass- weighted ro- 
tational velocity (L/r) versus enclosed gas mass 
(top) and the ratio of rotational velocity to cir- 
cular velocity versus enclosed gas mass (bottom) 
for the Zl— Gl and Z2— Gl runs. In both cases 
the typical rotational speed is 2—3 times lower 
than the circular velocity which means that the 
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Fig. 12. — Evolution of the spin parameter of the 
main halo for a metallicity of Z/Zq=10~^. 
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collapse of the halo is no t delayed by rotational 
support (|Wise et al.ll2QQ8f ). 
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Fig. 13. — Mass vs. radius plots for metallicities 
of Z/Z© = IQ-i (top) and Z/Z© = IQ-^ (bottom) 
for a radiation field of Go = 10~^ at redshift z = 
5. The lines represent the total (black solid), gas 
(green dashes), and dark matter (blue dot-dashes). 



3.4. Star Formation and Feedback 

When we include star formation and mechani- 
cal feedback in our simulations we see that, in the 
case of Z/Zq = 10 ~^ pre-enrichment, the halo is 
further enriched up to Z/Zq = 10~^ by the SNe of 
the first stars as early as 2: ~ 18. In Figure [HI we 
present the metal fraction (top) and the velocity 
divergence (bottom) of the central 10 kpc around 
the halo a short time after the first Pop III SN ex- 
ploded. In the velocity divergence plot, we see a 
clear ring-like structure with a large velocity gra- 
dient. This is due to the multiple outfiows driven 
by SNe that collide with infiowing material. Also, 
the highest outfiow velocities are right behind the 
shock front as expected. 

In Figure [161 where we plot time snapshots of 
temperature, density, and metallicity fraction of 
the central 25 kpc right after the first SNe, it is 
seen that the gas is heated to > 10^ K and expelled 
from the halo and hence infall of the outer halo gas 
onto the central core of the halo is quenched. The 
distribution of metals, provided by SNe, is very 
patchy. In the center of the halo, there are regions 
present with high metallicities, Z/Zq = 0.1, while 
in the outskirts of the halo the gas is still poor in 
metals. This indicates that SN feedback efficiently 
suppresses star formation in our enriched halos, 
irrespective of the ambient UV field. 

On the other hand, after the first SNe go off we 
do not see any significant difference in the distribu- 
tion of metals in the halo for the pre-enriched runs 
with Z/Zq = 10~^ and Z/Zq = 10~^, as shown in 
Figure [161 This indicates that the higher densities 
in the high-metallicity run do not have any effect 
on the mixing of metals in the halo. There might 
be a cooling instability effect on the mixing of the 
metals produced by the SNe with the surround- 
ings that we are not taking into account, but since 
the SN feedback is strong we think that this effect 
is negligible. 

4. CONCLUSIONS AND CAVEATS 



We summarize our results as follows. 
1) The redshift at which a multi-phase ISM is es- 
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Fig. 16. — From left to right: Slices of temperature, density, and metallicity fraction of the central 25 kpc for 
pre-enrichment metallicities of Z/Zq = 10~^ (left column) and Z/Zq = 10~^ (right column) for a radiation 
field of Go = 10~^ at redshifts z = 18.7, 15, 13.4, and 12.3 (from top to bottom, respectively) in the y — z 
plane. 
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Fig. 15. — Slices of metallicity fraction (top) and 
Fig. 14.— Mass-weighted rotational velocity velocity divergence (bottom) of the central 10 kpc 

(L/r) vs. enclosed gas mass (top) and the ra- ^ halo for Z/Z© = 10 at redshift 18 in the 

tio of the rotational velocity and circular veloc- U ~ ^ plane, 

ity against enclosed gas mass (bottom) plots for 
metallicities of Z/Zq = 10~^ (red dashes) and 
Z/Zq = 10~^ (black dots) at redshift z = 5. 
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tablished depends on metallicity. In the ab- 
sence of UV radiatio n, this critical metahicity 
is consistent with B romm fc Loebl (jloOS) and 
ISmith fcSigurdssonI ^20071 ). and Is {Z/ZQ)crit 



10 



-3.5 



2) Above a metahicity of 1% Solar the cooling 
efficiency of ambient gas no longer increases 
with a rise in metallicity (ISpaans fc Meijerink 

and therefore the Jeans masses in halos 
that are pre-enriched to metallicities of Z /Zq 
= 10~^ and Z/Zq = 10 ~^ are comparable. 

3) The cold dense gas phase is fragile to UV radi- 
ation for metallicities of Z/Zq < 10~^ and ro- 
bust to UV radiation for metallicities of Z/Zq 
> 10~^. Thus, the metal-poor star- forming 
ISM is fragile to UV radiation, and inclusion 
of a constant radiation background raises the 
critical metallicity value for the Pop III— Pop 
II transition from Zcr ~ 10~^'^Zq to Zcr ^ 
10-^Zq when Fq > 10"^ erg s ^ cm ^, which 
is in good agreement with the values tha t are 



found for the suppression of II2 by Shang et al 



(|2010h . This is because Ho is one of the main 
drivers of ion— molecule chemistry in metal- 
enriched gas and without it cooling is domi- 
nated by the fine-structure lines of [C II] and 
[O I] rather than CO. 

4) All the indicators for the dynamical evolution 
show that the halo evolves dynamically faster 
for Z/Zq > 10~^ and that the cold and dense 
gas phase in high metallicity halos survives (vi- 
olent) mergers due to the enhanced cooling. 

5) Pre-enrichment does not affect mixing of the 
metals that are produced by the first SNe. 

For future work, our simulations can be improved 
in several ways. We reach down to a maximum 
resolution of 75 and 7 pc with and without star for- 
mation feedback in our simulations, respectively. 
This is not enough to resolve the smallest gas frag- 
ments that will eventually turn into stars. Our 
simulations also assume a constant ambient UV 
background field which might not be realistic be- 
cause the UV radiation will originate from specific 
stellar sources. These sources are most likely to 
form in high-density areas and it would therefore 
be better to tie the UV background radiation to 



these sources at the time of their formation. 

Another improvement can be made by adjust- 
ing the way we pre-enrich the simulations. Cur- 
rently we assume one metallicity for a specific sim- 
ulation. We do this so that we can confidently 
study specifically the metallicity effects on the evo- 
lution of our halo. However, the distribution of 
metals, that are formed by the first SNe, is not uni- 
form but inhomogeneous. Therefore, in the near 
future we will perform simulations that consider 
patchy pre-enriched halos. 
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